Hanbury Brown and Twiss Correlations of Cooper Pairs in Helical Liquids 
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We propose a Hanbury Brown and Twiss (HBT) experiment of Cooper pairs on the edge channels 
of quantum spin Hall insulators. The helical edge channels provide a well defined beam of Cooper 
pairs and perfect Andreev reflections from superconductors. Surprisingly, the Cooper pairs do not 
bear a bosonic nature at all once emitted from the superconductor, and the cross correlation is 
strongly negative as for free electrons. This result is counter-intuitive, and seemingly contradictory 
to the previous theoretical and recent experimental results. We will point out that our setup is 
much closer in spirit to the original HBT experiment (and different) than those in previous works. 
We will argue that these differences affect significantly the results. 



The Bose-Einstein condensation of constituent parti- 
cles results in superfluidity in bosonic systems. 1 Super- 
conductivity in electronic (and hence fcrmionic) systems 
is attributed to the pairing of electrons into so called 
"Cooper pairs". 2 It is thus plausible to expect Cooper 
pairs to bear some bosonic nature, and to regard super- 
conductivity as a kind of condensation. Indeed, super- 
fluids and superconductors share many common proper- 
ties. Nevertheless, mathematically, Cooper pairs are not 
pure bosons because the pair creation and annihilation 
operators do not obey strictly the boson commutation 
relations, 2 and the bosonic nature of Cooper pairs still 
remains controversial. Therefore, it will be invaluable 
to test the issues directly in experiments. Here we pro- 
pose a Hanbury Brown and Twiss (HBT) experiment of 
Cooper pairs on the edge channels of quantum spin Hall 
insulators, 3, a recently discovered new state of matter. 
Surprisingly, the cross correlation is always negative and 
provides no hint of the bosonic nature of Cooper pairs 
once emitted from the superconductor. This result is 
counter-intuitive, and seemingly contradictory to the pre- 
vious theoretical 5-8 and recent experimental 9,10 results. 
We will point out that our setup is much closer in spirit to 
the original HBT experiment (and different) than those 
in previous works. We will argue that these differences af- 
fect significantly the results. Interestingly, the cross cor- 
relation has been predicted always negative in a diffusive 
multi-terminal superconductor-normal-metal contacts. 11 

A HBT effect is an intensity interference between two 
partial beams. It was originally introduced in order to 
overcome the technical difficulties in measuring the size of 
stars with Michelson interferometers. Since the pioneer- 
ing experiment by HBT in 1956, 12,13 it was soon realized 
that the effect can determine quantum-statistical prop- 
erties of a stream of particles: 14,15 The intensity correla- 
tion is positive for bosons (obeying Bose-Einstein statis- 
tics) while negative for fermions (obeying Fermi-Dirac 
statistics). 16 

An ideal HBT experiment requires a well-defined beam 
of particles and a tunable beam splitter. For normal elec- 
trons, an electron beam is realized by the chiral edge 
states of an integer quantum Hall insulator (IQHI), for 



which a quantum point contact (QPC) serves as a beam 
splitter. 17,18 It does not work for Cooper pairs because 
the high magnetic field required for IQHI breaks them. 

Quantum spin Hall insulator (QSHI) is a prototype 
example of a wider class of exotic states of matter, topo- 
logical insulators. Such states are characterized by bulk 
excitation energy gap and gapless edge modes, the latter 
being intimately connected to the topological nature. 3,19 
In QSHIs, the spin-orbit coupling gives the edge modes a 
helical feature, spin up (f) electrons moving to the right 
and spin down (J,) electrons moving to the left, and the 
time-reversal invariance prevents back scattering of these 
modes even in the presence of disorder. 3,19 The helical 
edge states of a QSHI are therefore a Kramers pair and 
can be regarded as a double copy of a chiral edge state. 
Then it is clear that a QPC serves as a beam splitter for 
helical edge modes. 20,21 

When a superconductor are put on top of helical edge 
modes, the edge modes become superconducting due to 
the proximity effect. 22,23 The edge modes beneath the 
superconductor will be referred to superconducting edge 
modes or simply "superconductors" . Then there forms 
a (point) interface between normal and superconduct- 
ing edge modes. A Cooper pair, itself being a pair- 
ing between time-reversal counterparts, moves with no 
back scattering across the interface between supercon- 
ductor and helical edge modes. That is, the Andreev 
reflection is always perfect even in the presence of inter- 
face barrier. 24 This is in sharp contrast to usual normal- 
superconductor hybrid structures, where in reality nor- 
mal electron reflection is unavoidable. Therefore, the he- 
lical edge states provide a well-defined and transparent 
channel for Cooper pairs. Two-dimensional QSHI has 
been observed recently in HgTe/CdTe quantum wells 4 
and also expected in InAs/GaSb/AlSb type-II quantum 
wells. 25 

Several HBT- type experiments of Cooper pairs have 
been proposed before. However, their setups all in- 
cluded a y-junction, instead of a true electron beam 
splitter, and not only Andreev but also normal-electron 
reflections are involved in the cross correlation. 7 It will 
be shown below (and in the supplementary) that these 
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FIG. 1. The schematic setups for HBT experiment of Cooper 
pairs using edge states of QSHI. The light gray represent the 
QSHI sample, the arrowed thick lines represent the helical 
edge states at its boundary, f (\.) electrons move counter- 
clockwise (clockwise). The QPC form by constriction with 
side gates serves as an electron beam splitter. A HBT ex- 
periment setup (a) for normal electrons and (b) for Cooper 
pairs. In (c), Cooper pairs are injected from two supercon- 
ductors. In this case, there occur multiple reflections between 
superconductors and the QPC. 
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FIG. 2. (a) Andreev reflection at the interface between normal 
(left) and superconducting helical edge modes (right), where 
t electron is reflected as a i hole, (b) The energy- momentum 
relation of the quasi-particle solutions to the Bogoliubov-de 
Gennes equations. Here c a (a =7, \) denote electrons and 
b a holes. The arrowed semicircles denote possible Andreev 
reflections. 



differences affect significantly the results. Our setup is 
much closer in spirit to the original HBT experiment. 
Very recently, the splitting of Cooper pairs has been 
demonstrated experimentally. 9 ' 10 However, in these ex- 
periments the scattering of Cooper pairs themselves have 
been suppressed with Coulomb interaction. While the 
experiments are only conductance measurement as they 
are, the cross correlation is therefore expected manifestly 
positive in such setups. 26 

We consider a stripe of QSHI. The helical edge states 
at its boundary are described by a one-dimensional Dirac 



H = hv 



E 



dx 



(1) 

where V't (tPl) is the field operator of "J" (4-) electrons, \x 
is the chemical potential, and v is the propagation ve- 
locity. The edge of a QSHI sample forms a closed loop 
and can never be terminated in the middle. In Eq. (1) 
it has been assumed that the lower (1=1) and upper 
(£ = 2) segment of the whole edge are isolated from the 
rest by putting contact reservoirs A, B, C, D, and other 
additional contacts, as shown in Fig. 1 (a). 

Constricting the QSHI bar with side gates forms a 
QPC. The f electron injected to the QPC from reser- 
voir A, for example, either moves on to reservoir C with 
probability amplitude t or tunnels to the upper edge go- 
ing out to reservoir D with probability amplitude r. The 
QPC thus serves as a tunable beam splitter for electrons. 
The QPC is thus characterized completely by the scat- 
tering matrix E 



(2) 



relating the amplitudes c' t(J (£ = 1,2 and a =t,|) at the 
output ports to ci a at the input ports of QPC. We have 
suppressed the energy (or momentum) dependence of q ct 
for notational simplicity. 

We apply a bias voltage V to reservoir C, keeping A, 
B, and D electrically grounded. With this bias config- 
uration, only 4. electrons are injected from C along the 
lower edge. The quantum statistical properties of the 
electron beam are characterized by the correlation func- 
tions defined by 
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dt (AI a (t)AI p (0) + AI (Q)AI a (t)) (3) 



with A I a (t) = I a (t) — I a , where I a is the average current 
into reservoir a. Given the scattering matrix in (2), the 
calculation of I a and S a p is a simple application of the 
Landauer-Biittikcr formalism. 17 We focus on the zero- 
temperature limit, ksT <C eV. In the present case, the 
average currents are given by 



e 2 V 
2nh 



e 2 V 
2tt^' 



(4) 



Clearly the injected current I = —Iq flows either to A or 
B (Ia+Ib — T)t but not to D (To = 0). This partitioning 
at the beam splitter leads to the well-known results 17 



SUa = S B b = -Sab = 2el\rt[ 



(5) 



for the current correlations, and its negative value is well 
understood by their antibunching behavior. 17 For later 
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reference, we stress that the contact D is completely 
"mute", i.e., Id = and Sad — Sbd = Sdd = 0. 

Let us now turn to the setup, Fig. 1 (b), of our main 
concern. We replace the normal contact C by a super- 
conducting contact S. A superconducting edge state is 
achieved by putting an s-wave superconductor on top of 
the edge. Due to the proximity effect, the edge states be- 
neath the superconductor form a superconducting state 
with an induced gap A = Aoe -1 ^. 23 The Hamiltonian of 
the superconducting edge modes is then given by [cf. (1)] 



Hs = hv I dx 



(6) 



We apply bias voltage V to S and keep A, B, and D 
grounded, effectively injecting Cooper pairs from S along 
the lower edge channels. The injection of Cooper pairs 
can be equivalcntly described by the Andreev reflection 
process, where say an f electron injected towards super- 
conductor is reflected by a 4- hole away from superconduc- 
tor, see Fig. 2. The Andreev reflection amplitude is given 
by a = e %{p ~ % arccos ( s / A o) f or incident electrons of energy 
E [\E\ < Ao). Note that due to the time-reversal invari- 
ance the Andreev reflection probability remains unity, 
| a. | 2 = 1, regardless of the imperfections at the edge- 
superconductor interface. This is in stark contrast to the 
usual normal-superconductor interface, where interface 
imperfections suppress Andreev reflections. 

Combining the Andreev reflections at the normal- 
superconductor interface and the normal scattering (2) 
at the QPC, one obtains the total scattering matrix for 
electrons and holes: 
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where bp„ and b 



(7) 

=t, I) arc amplitudes 



= 1,2 and 

of the holes on the input and output ports. Here r* = 
[r(-E)]*, t* = [t(-E)]*, and a* = [a(-E)]* describe the 
scattering of holes. We ignore weak energy dependence 
of r and t. 

The average currents I a and the current correla- 
tions S a p (a, (3 — A,B, D) are calculated using the 
Landauer-Buttiker formalism extended to the normal- 
superconductor hybrid structure. 2 '' 28 The currents 



T 2&2V \t\ 2 



In = 



2e 2 V 
2nh ' 



(8) 



are twice larger than (4), demonstrating perfect Andreev 
reflections at (or injection of Cooper pairs from) S. The 
correlation functions are also given exactly in the same 
form as (5), 



except that the total current I = Ia + Ib is now twice 
larger. Surprisingly and more importantly, the cross cor- 
relation is negative. (As in the normal case above, the 
contact D is mute; Id = and Sad = Sbd = Sdd = 0). 

Why is it surprising? The description in (7) in terms of 
electrons and holes is equivalent to Cooper pairs injected 
from S and scattered at QPC. Note that either (i) entire 
Cooper pairs go to A or B [Fig. 3 (a) and (c)], or (ii) 
constituent electrons in each pair split up into A and B 
[Fig. 3 (b) and (d)]. Naively, one may expect a positive 
contribution from case (ii) with one electron at each port 
A and B simultaneously. Assuming (partial) bosonic na- 
ture of Cooper pairs, one may also expect Sad > in 
case (i). 

The above naive expectation fails because it has ig- 
nored the two-particle interference. 29 For example, the 
two process (^4 f-electron; B ^-electron) — > (^4 |-holc; 
B |-hole) and (A f-electron; B f-electron) — !• (B |-holc; 
A 4_-hole) are not distinguishable and interfere with each 
other. Due to the Fermi-Dirac statistics, the amplitudes 
for these processes are opposite in sign, hence giving neg- 
ative cross correlation in (9). 

There is another simple way to understand the nega- 
tive correlation. For example, f electron from A under- 
goes either Fig. 3 (a) or (b). Since the perfect Andreev 
reflection at the edge-superconductor interface is noise- 
less, it is nothing but the partitioning of a single hole and 
thus gives a negative contribution to Sab ■ Similarly, the 
processes in Fig. 3 (c) and (d) also give negative contri- 
butions. 

We thus conclude that Cooper pairs bear no bosonic 
characters at all, once they get out of the superconductor. 
It is interesting to note that in recent experiments with 
quantum dots 9 ' 10 the processes Fig. 3 (a) and (c) have 
been suppressed due to Coulomb interactions. In this 
case, there is no two-particle interference. 

The above conclusion appears contradictory to the pre- 
vious theoretical works, 5-8 where positive correlation was 
predicted in certain range. As mentioned earlier, these 
works all include a Y junction, not a true electron beam 
splitter. One crucial effect is the multiple reflections be- 
tween the junction and the superconductor. Further, not 
only Andreev but also normal-electron reflections are also 
involved. The latter effect cannot be simulated in our 
system and is analyzed in the supplementary. 

To simulate the former effects, let us now replace both 
C and B by superconductors, <Si and S2 , with phases </?i 
and ifi2, respectively. We apply bias voltage V on both 
Si and S2 and keep A and D grounded. The scattering of 
electrons and holes is governed by the scattering matrix 
(ignoring corrections of order E 2 / Aq -C 1) 
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Saa = S BB = Sab = 2el\rt\ 2 , (9) 
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FIG. 3. Elementary scattering processes in a Hanbury Brown 
and Twiss experiment of Cooper pairs, described in electron- 
hole picture. Blue (red) line represents the edge channels 
of f (i) electrons (filled circle) or holes (empty circle). In 
Cooper pair picture, panel (a) corresponds to Cooper pairs 
being injected from S, scattered at the quantum point contact 
(QPC), and going entirely to A. Analogously, in (c), Cooper 
pairs entirely go to B. In (c) and (d), constituent electrons 
in a Cooper pair split up at the QPC. 

where a\ = -e !f,1+ ' af-il , a 2 = _ e 'V2-i2*f i 2) kp 
// /;</. and L\ and Li are the distances from the QPC 
to the superconducting contacts S\ and S2, respectively. 
The multiple reflections between QPC and superconduc- 
tors are manifested through the factors z = 1 — \r\ 2 e te 
and w = 1 — e l6 , where 9 = (<pi — y> 2 ) + 2kp{Li + L%) 
is the phase accumulation during one cycle of the multi- 
ple reflections. This immediately leads to the resonance 
behavior of the average currents, 

- - 2e 2 V \t\ 4 

Ia ~ Id - ^Tl + |r| 4 -2|r| 2 cos0 (U) 



and the correlations 

Saa = Sun = Sad = 2el- - ^ cos9 (12) 

For 8 = 2im with n an integer, the currents take maxi- 
mum and the noises vanish. 

We note that the sign of the cross correlation Sad 1S 
now positive, in contrast to (9), which is always nega- 
tive. This difference is ascribed to the multiple reflec- 
tions between QPC and superconductors. As noted by 
Baym, 30 the intensity correlation (I(t + r) I(t)) is equiv- 
alent to the relative probability to observe two particles 
at two points separated by a distance vt, where v is the 
flight velocity. If the HBT experiment is done with a 
"true beam splitter" as in Fig. 1 (b), then the intensity 
correlation is entirely due to the spatial distribution and 
the quantum statistical property of the particles in the 
beam itself. To the contrary, if the beam splitter is re- 
placed by a Y-junction and multiple reflections occur 
between the source and the junction, then the cross cor- 
relation is not directly related to the spatial distribution 
in the beam alone but affected significantly by the succes- 
sive interactions with the source. Therefore, the positive 
cross correlations in Refs. 5 ~ 8 do not represent entirely 
the quantum statistical properties of Cooper pairs after 
emitted from superconductor (away by distance larger 
than the superconducting coherence length). 
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